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Abstract 

In a recent paper, Creutz has given a new action describing two species of Dirac 
fermions with exact chiral symmetry on the lattice. This action depends on param- 
eters which may be fixed at certain values in order to get the right continuum limit. 
In this letter, we elaborate more on this idea and present an action which is free of 
any other parameter except the fermion mass. 

PACS numbers: ll.15.Ha, 11.30.Rd, 71.10.Fd, 71.20.-b 
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Ever since the birth of lattice QCD it was understood that chiral symmetry cannot 
be easily realised on the lattice pQ. The naive discretization of the Dirac operator on a 
regular hypercubic lattice, as given in the momentum space, 



possesses exactly 16 zeros in the 16 corners of the Brillouin zone. Since they differ only 
on their slope, one can interpret them as 16 degenerated Dirac species. 

Later, Wilson lifted this degeneracy by defining a new operator, namely [2], 



This way, one gets a single Dirac fermion on the lattice with the chiral symmetry explicitly 
broken. 

Another approach, the staggered fermion approach, puts the four components of a 
Dirac spinor on different lattice sites [3]. This reduces the degeneracy to four species. 

Further research on chiral symmetric lattice fermions was discouraged by Nielsen and 
Ninomiya theorem, which states the impossibility of a single left handed chiral fermion on 
the lattice [I]. However, as pointed out by Wilczek, the theorem does not prevent having 
a pair of Dirac fermions on the lattice [5]. So, he argues, if flip), f2(p), fsip), fiip) are 
suitable functions of lattice momenta one can get a Dirac operator, 





with only one doubler. His explicit construction, 




A > 1 



1 

- sin PjQ, 



J = 1,2,3, 



2 



exhibits two zeros: one at (0,0,0,0) and the other at (0, 0,0, 7r). However, the loss of 
cubic symmetry introduces extra relevant terms in the Wilson plaquette action [5]. 

Recently, motivated by the Dirac structure of graphene electrons in two dimensions, 
Creutz was able to elegantly generalise this structure to four dimensions, exactly what 
one needs in particle physics [6]. This generalisation inherits the gapless property of the 
electronic pair which comes from the fact that the electron and its doubler sit symmetri- 
cally with respect to the centre of the Brillouin zone. He starts by defining the left handed 
degrees of freedom in 4-momentum space: 



where er^, k = 1, 2, 3 are the Pauli matrices and B, C are real parameters. In order to get 
the Dirac fermion we bring this expression in the form: 



zip) = B [AC — cospi — cosp 2 — cosp 3 — cosp 4 



+za 1 (sinp 1 + sinp 2 — sinp 3 — sinp 4 ) 



+za 2 (sinpi + sinp 2 — sinp 3 — sinp 4 ) 



+icr 3 (sin pi + sinp 2 — sinp 3 — sinp 4 )] 
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where Skip) are defined through: 



s iip) = srn Pi + sinp 2 — sinp 3 — sinp 4 



s 2 (p) = sinpi — sinp 2 — sinp 3 + sinp 4 • 



Sz(p) = smp 1 — smp 2 + sinp 3 — sinp 4 



Introducing right movers as z(p) we get the Dirac operator in the form: 



/ 



\ 



zip) 



Zip) 
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Using the Euclidean Dirac gamma matrices in the chiral representation 



/ \ 

-ia k 



Ik 



( 



74 



/ 



we get the momentum space Creutz-Dirac operator: 



aD c {p) = iB^ A (AC - Y cosa Pp.) + iB ^J7 fc s fc (ap) 



k=l 



where we have restored the lattice spacing a. Creutz has shown that the zeros of this 
operator in the Brillouin zone are at (p,p,p,p) and (—p,—p,—p,—p), where C = cos p. 
The reader can immediately see that the formal continuum limit, i.e. a — > 0, does not 
give the continuum Dirac operator j>. While this is not a problem, it is related to the 
fact that the Creutz action has no zero at the origin, making its continuum limit subtle. 
In order to get a hyper-cubic lattice of p-translations, Creutz shows that its parameters 
must satisfy BS = C, where S = sin p. 

From the discussion above it is clear that it is desirable to have an action with or- 
thogonal p-axes which has a formal continuum limit. While these features may not be 
necessary [7j, they serve us as a guiding principle to fix the parameter C. We shall show 
that such an action exists and has the form: 



d (p) = ^ sin *v + Yl i %( cos - x ) 



1.21 



where 7' \x = 1, 2, 3, 4 is another set of Dirac gamma matrices, which are given as a linear 
combination of 7^, 

V 

and the coefficients of the orthogonal matrix a. This action possesses one zero at 

p = 0, the location of the rest depending on the details of 7'-matrices. In the rest of the 



paper we will show that this action can be derived from the Creutz action. In doing so 
we will find the concrete form of a. 

We start by dropping the overall parameter B and setting again lattice spacing to one: 



D c (p) = ?74(4C cos ap^) + % ^ikSkM . 

\i k=l 

Translating momenta, p^ = p + q^, we have: 

3 3 

D c (q) = ?74(4C - C ^ cos %) + iC ^2 ^kSk(q) + ^74 ^2 sin % + iS ^2 ^Ck{q) ' 

H k=l pi k=l 

where 

ci(q) = (cosgi - 1) + (cosg 2 - 1) - (cosg 3 - 1) - (cosg 4 - 1) 
°2(q) = (cosgi - 1) - (cosg 2 - 1) - (cosg 3 - 1) + (cosg 4 - 1) • 
c 3 (q) = (cosgi - 1) - (cosg 2 - 1) + (cosg 3 - 1) - (cosg 4 - 1) 
In order to simplify the action we assume C + S = 0, which gives C = l/y/2 conforming 
the C > 1/2 constraint imposed by Creutz. Rescaling the action and defining: 

s i(q) = — sin g x — sin g 2 — sin g 3 — sin g 4 

c 4 (g) = (cosgi - 1) + (cosg 2 - 1) + (cosg 3 - 1) + (cosg 4 - 1) , 

the translated Creutz action takes the form: 

D c (q) = ^ilnS^q) + ^2h^c^q) , 
n n 

or in the scalar product notation, (7, x) = X^T/i^M' one has: 

D C (q) =i(l,s(q) + c(q)) . 



5 



Now, introducing the following orthogonal matrices: 



1 



1-1-11 
1-11-1 



v -i -i -i -i, 



6: = 



1-1-1 1 
1-11-1 

v 1 1 1 V 



and noting that, 



s = 2as, c = 26c , 



where 



s = ( sin^i, sing 2 , sing 3 , sing 4 ) T , c=(cosgi — 1, cosg 2 ~l, cosg 3 — 1, cosg 4 — 1 ) T 



then, the rescaled action by a factor of 2 can be written in the form: 
D c (q) ■= 1(7, aS(g) + 65(g)) = i(a T 7, s(q) + a T bc{q)) 

Denoting, 

/ 

a := a b 



-1 
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1 


V 
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-v 



we get: 

D c (q) = «(a T 7, s(p) + atc(q)) . 
It is easy to show that a T 7 are Dirac gamma matrices □ , therefore the factor a T can be 



droped. This way, our final expression is: 

Dc{q) = i{l,s(q))+i{i,c{q)) 

= Yl ^ sin + 2 Z7 m( cos Qi* ~ x ) ' 

where 7' = 07 are again Dirac gamma matrices. This is exactly the action of eq. (jl.2p . 
Noting that, 

A 4 A* 

we get another expression for the fermion action: 

D c(q) = £ i'7/i sin 9m + £ ^cos g M - 2C . (1.3) 

n n 

It is easy to see that the formal continuum limit of this action is: 

Dc{q) -> ^ , 

hence, this action has one zero at q — 0. It is easily verified that the other zero is at 
(7r/2, 7r/2, 7r/2, 7r/2). Since we arrived at this action starting from the Creutz action with 
equivalent transformations, which conserve the number of zeros, then we may conclude 
that there are no more zeros. Hence, the action defined in eq. (jT3j) describes two species 
of chirally symmetric Dirac fermions. The slope at the second zero can be computed by 
expressing g M = n/2 + in eq. (II. 3ft : 

D c (tt/2 + k) = H» (cos fc M - 1) - £ ^sin k^ -> - ^ i^A; M = - ^ 7 M (aA;) M . 

A* A 1 A* A* 

Therefore, the slope at (tt/2, tt/2, tt/2, 7r/2) is the matrix —a. 

In order to write down the action in the position space, we use eq. (II. 3p . Adding a 
mass term and multiplying both sides by i we get: 

iD c (q) = im + 2Y- l -Y,[^- ^) ^ + (7; + *'7j e^] • 

A* 
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Hence, the position space Dirac operator involves the following three terms: 
the diagonal term, im + 21"; 

the forward hopping term in the /i direction, — |(7u — 

the backward hopping term in the \i direction, — 1(7„ + &7,u)- 

Finally, the gauge fields can be introduced in the standard way, the resulting action being 
gauge invariant. 

In conclusion, based on the original action of Creutz, we have given a minimally 
doubled and chirally symmetric fermion action. The formal continuum limit served us 
as a guiding principle to fix its parameters. However, the parameters have geometric 
significance and can be fixed at other values in order to maximize the symmetries of the 
action 0,[7]. As was shown in detail at references [3 IE], the lack of hyper-cubic symmetry 
may introduce relevant terms in the presence of interactions. Thus, as we already know, 
chiral symmetry on the lattice comes with a high price. Whether the price of the minimally 
doubled action is higher or lower than the price of other actions, this is a question that 
can be settled by direct simulations [9]. 
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